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Abstract 



In this paper we provide a detailed numerical study of the influence of 
thermal effects on the original picture of the Affleck-Dine (AD) baryogenesis. These 
effects are found to modify the results greatly in some cases. We estimate the 
• baryon/entropy ratio and provide numerical results on the typical behaviour of the 

charge as a function of the strength of the potential and other parameters. 
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1. Introduction 



AfHeck-Dine baryogenesis from flat directions [1, 2] is a natural mechanism to ex- 
plain the baryon/entropy ratio. The current known value of Q^/i^ is 0.019 ± 0.002 from 

nucleosynthesis and 0.031 ± 0.005 from BOOMERANG data [3]. This corresponds to 
ns/n^ ~ 5.1 X 10^^ and nB/n^ ~ 8.3 x 10^^ respectively. The ingredients needed in 
general to generate the baryon asymmetry are CP violating A-terms, terms which lift 
flat directions, and supersymmetry breaking terms in the early Universe which give rise 
to negative mass squared of order of if^, where H is the Hubble constant. At the level of 
renormalizable interactions there are several flat directions which can be found in MSSM, 
i.e., directions where F and D terms in the potential vanish. A simple example is the 
HuL flat direction. These flat directions can be lifted by nonrenormalizable terms which 
we might imagine to be associated with the Planck scale, i.e., suppressed by some power 
of M*. ^ In that situation as a result of the balance between soft induced SUSY breaking 
potential and higher order nonrenormalizable terms the field can acquire a large expec- 
tation value which evolves as some power of t. Eventually, the negative mass squared 
term, which decreases as 1/t^, becomes comparable with 777.3/2 term and the AD- field 
starts to oscillate. This is the moment when baryon charge, which is generated by the 
torque due to the different phases and time dependence of A-terms, "freezes" at some 
value. Terms which lift the fiat direction in general have the form ^frr- We will prin- 
cipally consider 77 = 1, 2, 3 and refer to these as the n = 1, 2, 3 cases. One then can do 
a simple estimate to obtain the ratio ns/n^^ . For example, in the n = 1 case, ub/ti../ is 
of order (7773/2/M*)^/-^sin((5) ~ 10^'^sin((5) or so. ^ The resulting baryon/entropy ratio is 
a bit lower because one has to take into account that, for example, in the case of the 
H^L direction lepton number is produced first, which later has to be converted to baryon 
number. Also, the AD condensate interacts with thermalized infiaton decay products 
and some part of it can be evaporated before the charge is produced. There is a related 
issue of Q-ball formation and the corresponding evaporation rate, which is different from 
that of homogeneous condensate. There have been a number of papers on the subject in 
which authors considered thermal effects relative to AD baryogenesis as well as Q-ball 
formation. In these paper we will consider AD baryogenesis in the absence of Q-balls. 
Taking into account the formation of Q-balls will alter the whole picture. We refer the 
reader to [4] as for this subject. 

Recently it was noticed [5] that due to the evolution of the AD condensate in the 
background of infiaton decay products there is an interesting effect that takes place in 
addition to the physics of the original scenario. It was observed that superpotential in- 
teractions couple the fiat directions to other fields. These fields acquire masses induced 
by the flat-direction vev but they may be sufliciently small so that flelds come to thermal 

-'^By M* wc denote reduced Planck mass, M* = lO-'^^GeV. 
^ By (5 we denote CP- violating phase. 
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equilibrium with the inflaton decay products. In such cases the flat direction starts to 
oscillate at earlier time than usually estimated because it acquires thermal mass y'^T'^^ 
which decreases with time as 1/t. Since —H"^ falls more rapidly with time, the differ- 
ence y^T^ — H"^ eventually becomes positive. That normally happens much earher than 
t ~ 7^3/2. 

It was also argued that the main source for generating baryon asymmetry in that 
case are A-terms which are proportional to T. In [6] it was pointed out that such A-tcrms 
are suppressed, in general, by symmetries so that there is no temperature enhancement of 
the A-teim. However, it was noted that there is an additional source of A-terms that can 
be efficient in the n — 2,3 cases. It was also shown that there is an additional thermal 
contribution to the potential of the form T^log(|0p). In the n = 2, 3 cases, this defines 
the time when the condensate starts to oscillate, rather than the y'^T'^ contribution found 
in [5]. 

In this paper we further investigate this scenario. We consider different types of 
thermal effects which are relevant for different choices of n. We analyse the parameter 
space in greater detail, i.e., the dependence of the resulting baryon/entropy ratio on the 
parameters of the lagrangian. These parameters are the Yukawa constant, gauge coupling 
constant, the relative phase between A-terms, and the coefficients in front of the A-terms. 

The n = 1 case is found not to generate sufficiently large nB/n.y, but the n = 2,3 
cases can rather easily generate the needed number for a wide range of parameters. 
Throughout wc assume for simplicity that the ratio of the inflaton mass to the reduced 
Planck mass is of order of 10^^, taking, in general, the mass of inflaton mj to be lO^^GcV 
and reduced Planck mass to be lO-'^^^GcV. These parameters appear in the estimate of 
the reheating temperature Tr as well as in the estimate of the baryon number violating 
terms. Therefore, the choice of these parameters is important for the estimate of ns/n-y. 
We discuss this in detail in Section 4. 

This paper is organized as follows: in Section 2 we discuss the origin of two thermal 
effects as well as their relevance in the n — 1,2,3 cases. In Section 3 we introduce 
additional A-terms which will be the sources for generating the baryon asymmetry. In 
Section 4 wc provide the estimates and numerical rcsTilts for the baryon/entropy ratio. 
In Section 5 we give some details of the numerical simulations. 



is the coupling constant and T is the temperature. 
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2. Thermal effects 



As was argued in [5] and [6] there are two types of thermal-induced contributions 
to the potential. The first is due to the following mechanism. Consider some field 
X which couples to the flat direction via the superpotential W ~ yxx4' interacts 
with the dilute plasma produced by the infiaton decay. Because of x's coupling to the 
flat direction, chi acquires mass m^^ = y4>. If this mass is less than the temperature 
of the thermal plasma this fleld will come to thermal equilibrium, giving an effective 
temperature-dependent mass to the flat direction. Let us to do some estimates to see 
when this effect might be important. Consider the case n — 1. Since T — (i/Ff^M^)^ 
and (p — {HM^y/"^, the condition yep < T leads to 

where Hth is the value of Hubble constant when thermalization of x occurs. Taking M* 
to be lO^^CeV and Tr = 10^°GeV one obtains: 

Htn < (^)^10i°Gev. (2) 

y 

On the other hand, if gets a thermal mass, oscillations start when yT ~ H, which gives 
the value of the Hubble constant Ho at this moment 

Ho ~ y'/'T'J'Ml/' = (^)^/^10^°Gev. (3) 

In order for this thermal term to be relevant oscillations have to start not later than x 
becomes thermalized, so that we have H„/Hth < 1, which means y < 0.01. Provided y 
is sufficiently small, there is a y^T^|<^p-contribution to the potential, which affects the 
time when oscillation of begins. If y is bigger than 0.01, then oscillations start when 
t/^T^ term is generated, which means that Ho is defined by Hth- Estimating Hth for the 
n — 2 case gives: 

n^2:Hth< ^ = (^)^^10-«GeV, (4) 

which is too late for thermalization process to have any effect on V^cj)) unless y < 0.001. 
For n — 3 from the condition ycj) < T one obtains (q^) < 10~^. We conclude that in the 
n — 2,3 cases this effect does not take place unless y is very small. 



The other contribution to the potential comes from the modification of the coupling 
constant when some sfermions which are coupled to the fiat direction gain masses. The 
effective potential for (see [6]) is then 

Veff ~ «^rMog(|0n, (5) 

^Frf is the inflaton decay rate. 
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where ol is the gauge couphng. This contribution causes to oscillate when ~ 
g^j/ ■ For n = 1, starts to oscillate at ~ giTr = (qw)-'-'^^^^^' ^'^^ n = 2: 
Ho ~ a6/5TK(^)V5 = (_ii_)6/5io6GeV, while for n = 3: ' i/^ ~ a4/3y^(^)V3 = 
(Q^)'^/^10^GeV. If we now look at the ratio of HgscS due to both thermal contribu- 
tions in the n — 1 case, we see that 

||.l^|..^103, (6) 

H^^> a Tr a 

where by Hj^^^ and if^^^ we denote the values of the Hubble constant when the AD 
field starts to oscillate if there is only Q;^T^log(|0p) or y^T"^ thermal contribution to the 
potential, respectively. 



That estimate shows that in the n = 1 case the y^T^ term dominates over the 
a^T^ log(|0p)-potential unless the ratio < 10"'^. Therefore, we will analyze the n = 1 
case with the y^T^ term only, while for n = 2, 3 we will consider only the logarithmic 
thermal term. 



3. A-terms 

To build up a baryon or lepton asymmetry one needs to have corresponding U{1) 
violating 74-terms. The misalignment between their phases then exerts a torque, making 
(f) to roll down towards one of the discrete minima along the angular direction and settle 
there. In the original scenario [2] there were two A-terms : 74777.3/20""*"^ /M", which is the 
usual supersymmctry breaking A-term, and ai/0"+^/Af", which is induced because of 
the finite energy density in the Universe during inflation.^ At later times {H < 171^/2), 
because of its dependence on H, the second term is no longer important, but at the times 
H ~ m3/2, when the flat direction starts to oscillate, this term is of comparable size with 
the first one and a sufficient amount of charge is produced. 



As we learned in the previous section, 0-oscillations generically start much earlier 
than H ~ 777.3/2. One might then consider some other sources of A-terms, which can be 
relevant in this case. These additional A-terms may arise in the following way: consider 
the superpotential 

sw^jd^em^, (7) 

where /(/) is some holomorphic function of ithe infiaton field /. Taking the first two 
and a are some complex constants. 
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terms in the polynomial expansion of /(/) in powers of I /M^ leads to the superpotential 



1 P 

5W = —{aI + b—)^, (8) 

where, in general, a and b are complex constants which do not need to have the same 
phase. / decreases as Io{to/t), so that the second term is somewhat suppressed with 
respect to the first, but at the time when oscillations of the (^-Held start this is not 
necessarily a large suppression and one might expect to generate a reasonable baryon 
number. We will investigate this particular case in greater detail in the next section. 



4. n^/n^: estimates and numerical results 

Before we move to the discussion of numerical results let us provide some crude 
estimates for the baryon number. In the n = 1 case we assume that y < 0.01 and y'^T'^ 
dominates over the logarithmic term. Then, Hg ~ 7/^/^Tj^(^)^/^. We take the equation 

^-\Vb\MS), (9) 

where 5 is some effective phase which comes from the phase difference of two 74-terms, 
Vb is baryon number violating part of the potential. For the A-tcrms which come from 
superpotential 5W = -^{al + bjj^)^^ S would be the phase of b after we redefine the 
phase of 4> to make the al term real.^ 

Taking the initial amplitude of / to be of order of M* and replacing by the 
product ubHo one obtains 

We compute n^t^ instead of because it is more convenient, since after the oscillations 
begin the baryon density decreases as 1/t'^. For the same reason it is more convenient 
to compute n^t'^ — T^t^ where ta — (Fd)""^, and F^ f« my Ml is the inflaton decay rate. 
Since T = (HTdMlY^^, the reheating temperature is given by Tr = (F^M^*)^/^. Then, 
T|t2 = Ml'^t]!'^. Taking mj ^ lO^^GeV one obtains F^ ^ 10''^GcV, Tr ^ lO^^GcV and 
n^t^ ~ lO^^GeV. The value of Tr is somewhat large from the perspective of the gravitino 
problem. The baryon/entropy ratio is 

n=l:^«6(^)'/'lO-^'sin(5,), (11) 



^The phase of / is assumed to be a constant, so that one can absorb its phase in b and consider the 
inflaton to be real. 
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which turns out to be too small unless y ~ 1. Actual numerical study shows that for a 
wide range (see Fig. 1) of y the resulting baryon/entropy ratio is somewhat larger: 

— ^ (bsm(Sb)){10-^* - 10-^^), (12) 

and almost independent of y. There is some uncertainty which depends on the choice 
of as well as m/. However, in order for n = 1 case to be viable we would need to 
consider a much smaller ratio of M*/m/ than 10^. 

One can also consider = emf/M^, where e is some new parameter, i.e., F^ = 
mf/A^, where A is some new scale different from M*. Then Tr is lower than before by 
a factor of e^/^. For e ~ 10^'^, ~ lO^GeV, which is an upper bound on Tr to avoid 
gravitino problem. If one combines the result of Eq. (10) and the value for T|t^ in terms 
of e, the inflaton mass, and M*, one gets that ris/n^ = hy^/^e'/^{mi/M^f'^. This IS an 
even smaller number than at Tr = lO^^GeV. So if one tries to fit the value of Tr in 
a consistent way to avoid gravitino problem one obtains that the n — 1 case seems to 
be even less acceptable. In the literature sometimes the value of Tr is considered as a 
free parameter. In that case the value of the ratio ^ ~ ''^^ is very sensitive to the 
value of the reheating temperature and our estimate for the resulting baryon/entropy 
ratio changes by several orders of magnitude. For example, if one takes Tr = lO'^GeV 
again and keeps F^ lO^GeV one will obtain ^ ~ 10~^. One can also treat td as a 
free parameter as well. This, in fact, would significantly relax constraints on nst^ and 
would soften our conclusions. There was also the discussion in the literature [7] , that the 
gravitino problem can actually be avoided even at temperatures higher than lO^^GcV. 
In this paper we stick to Tr ^ lO^^GeV according to the estimate which wc get from the 
expression T = (ifFj^M^)^/^ for the temperature during the inflation at if ~ F^ as well 
as the estimate for Fd mentioned above. Therefore, for that choice of parameters, we can 
conclude that n—1 case can hardly be acceptable. However, with all the remarks above, 
one can probably consider the n — 1 case to be a borderline case. 



In the case of n = 2, a similar analysis gives now there is no 

y^T^-term and oscillations start due to the logarithmic contribution). Since (j) behaves 
as [HMiy/^ one obtains the following estimate for net^: 

n = 2 : nBt" ^ ba^/'T^/'Ml/' « 6(^)^/^10^iGeV, (13) 

and the corresponding baryon/entropy ratio is 

n^2:'^ ^b(—Y/^10-'^sm{Sb). (14) 
0.1 

This is very different from the numerical result, which is 

— ^ (6sin((56))(10-^° - 10-^), (15) 
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for wide range of a (see Figs. 2, 3). We want to point out here that as in the case n = 1 
the dependence of the result on the couphng constant is rather weak, as opposed to what 
one gets estimating ns/n^ analytically. 

Before we explain why the numerical result is so different from the naive estimate, 
we want to repeat the analysis above for the n = 3 case. I that case (p = (HM^Y^*; 
the operator which creates the charge is ^^j^- Estimating n^i^ a,t H — Hg one gets 

net^ ~ bHy^Ml^^sm{Sb). Since ^ a'^l^'T^J^'/Ml'^ 

UBt" ^ ba''/^ sm{Sh)Tl/^Ml/\ (16) 
Taking Tr = 10^°GeV and = lO^^GeV one gets 

UBt^ ~ ba^/^ sm{Sb)10^^GeV, (17) 
and, the corresponding baryon/entropy ratio 

- = ^^ba'/'sm{5,)10-'^ (18) 

which is too small for all reasonable values of 6, y and sin((5b). One can notice that this 
rough estimate predicts the decrease of baryon/entropy ratio as y gets smaller as well as 
in the n — 2 case. Actual numerical study, as in case n = 2, gives a quite different value 
of that ratio. The baryon/entropy ratio for the wide range of a is 

— «6sin(56)(10-^-10-^), (19) 
which could be consistent with experimentally known value if 6 ~ (0.1 — 0.01) and ~ 1. 



The naive expectation fails because the approximation of dns/dt by nsHo is too 
crude with respect to actual integration. One can understand why the amplitude of n^t^ 
grows as y gets smaller, when looking at the behavior of rifit^ at some values of y as a 
function of time. Prom numerical results it is clear that the later the oscillations of 
start the more oscillations nst^ undergoes before it "freezes". With each oscillation, the 
operator which is responsible for the charge production contributes more and more to the 
amplitude of nst^ ■ That might explain the unexpected behavior of charge in the n = 2, 3 
cases. One can see that the behavior of the charge for n = 1 with the y^T^-term is quite 
different from its behavior for n — 2,3 with the logarithmic term. Namely, for n = 2, 3, 
before the value of a charge "freezes" , it experiences many more oscillations than in the 
n = 1 case, gaining larger amplitude with each of them. One can see this from numerical 
results in Fig. 1-6, which show the evolution of the charge with time and with the strength 
of the potential (Yukawa coupling in the case n — 1, and gauge coupling in the n = 2, 3 
cases). The first three plots demonstrate the expected behavior of the charge, i.e., it 
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oscillates at the time t < to and freezes afterwards. The difference between n = 1 and 
n = 2, 3 is that for n = 1 the charge typically oscillates one or two times, while in the 
n — 2,3 cases it oscillates much more. Correspondingly, for n = 1, the numerical result 
is close to the naive estimate, but for n = 2, 3 numerical value is much larger. 



The Figs. 1-3 further distinguish these cases. One can see that changing the Yukawa 
coupling in the phenomenologically interesting region does not produce major changes 
in the amplitude of n^t^ in the n = 1 case. However, if one changes the value of the 
gauge coupling in the n = 2, 3 cases one can sec that n^t^ oscillates with log(a) and its 
amplitude grows with the decreasing of a, which is against the naive analytical estimate. 
We believe that this happens due to the logarithmic nature of the thermal potential in 
these cases. 



Dependence of numerical result of the baryon/entropy ratio on b and 5b doesn't 
bring any surprises with respect to the crude estimate above. In Fig. 7 one can sec that 
dependence on 6i, is sin(5{,), and from Fig.8 the dependence on 6 is a linear function. Thus 
for any b and 5^ one knows how to reproduce the baryon/entropy ratio from the results 
in the equations (12, 15, 19). 



5. Details of the numerical simulations 

The equation for the evolution of the (f) is 

dV 

'/• + 3i^0+^ = O. (20) 

In general V has the form V = Vo(|0p) + Vb(0) + Vb(0*), where Vq determines the 
evolution of |0| and Vb is the part of the potential which is responsible for the creating of 
the baryon charge. It also has a minor effect on evolution of |0|, which can be neglected 
in the cases we consider. The potentials we wish to study have the form: 

Vo = i-H^ + m^/,)|<^p + V,,{T) + l^p, (21) 



Vb = aHW + b^W + Am^i<,W 

As we discussed before, the terms which are proportional to 777,3/2 are not important. 
Dropping them, for n = 1 we have: 

V = (-7/^ + y^r^m + ^% + (a/^^ + + 1..C.). (22) 

for n = 2: 

V = -E-\4>\- + a-T^ logd^r) + A^M! + («^,^ + ^'^^ + h.c.), (23) 
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and for n = 3: 



We took into account that different tfiermal effects are relevant for n = 1 and n — 2,3. 
Introducing radial and angular parts of by = Re^^ one gets corresponding equations 
for R and O: 
n = 1: 

R + 3HR + i-H^ - + ^2j.2 ^ A^)^ ^ 0, (25) 
n + {3H+ f )Q = |a|^ sin(4n) + sin(4n - 6^)] 



n = 2: 



n — 3: 



i? + 3HR + i-H^ - + Q,2T7i?2 + >^)R ^ 0, 
n + (3i/ + = |a|^ sin(5fi) + \b\^ sm{5n - 



(26) 



R + 3HR + i-H^ - + a^T^R^ + ^)R = 0, (27) 
h+{3H+ f )0 = \a\^sm{6n) + sin(61^ - S^). 

We are dropping the contribution from at t < t„ to the equations for R for the 
following reasons. Let us to compare the behavior of R induced by the and the 

|a|-^:^ cos(4f^) terms separately (n = 1 case). We know that the first term causes R to 
behave as {HM^Y^"^. The second term causes R to behave in the same way, provided 
VL is small, so that cos(4r2) pa 1. The smallness of Vt is guaranteed by its equation of 
evolution and, in fact, is well seen in numerical results. Therefore, qualitatively it does 
not bring any changes to the evolution of R7 The other term, which is proportional to 
is suppressed by (mj/M^,) and falls down more rapidly with time with respect to the 
A-term. Hence, neglecting these terms in the equation for R makes sense, since that does 
not change the evolution of R and simplifies the numerical and analytical analysis of the 
equations. 



For numerical simulations it is more convenient to introduce a dimensionless field 
r = Rt and scale everything by m/, which we take to be 10-^^GeV. Then, the equations 

^In fact, one can also neglect {l? term in the equation for R since this term is small compared to the 
others for any H > Ho- It becomes, however, important after the oscillations start. At this moment 
baryon charge is already produced and fixed. Later evolution of the baryon number is not interesting 
because of that and one does not need to run simulations further. 
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above with all the simplifications which were discussed before take the form: 
n = 1: 

r + {-H^ + y^T^ + e^^)r = 0, (28) 



Q + f Q = |a|e^ sin(4Q) + jfeje^^ sin(4Q - 5^); 



n = 2: 



n — 3: 



r + i-H^ + aH^ryr^ + e^^)r = 0, (29) 
Q + f n = |a|e'^ sin(5Q) + j^le^^ sin(5Q - Si,); 

r + {-H^ + aH^T^r^ + e6^)r = 0, (30) 
n + ffl^ |a|e3^ sin(60) + |6|e^^ sin(60 - 5^), 

where e = Mj/M^ 10~^. Solving these equations numericaly and noting that ns^^ = 
r^Cl, we can plot the behavior of n^t^ as a function of time as well as find its asymptotic 
value as a function of different parameters. 



At this point it is rather easy to analyze the behavior of the charge. In fact, the 
equation for r is no longer dependent on Q, so that the second equation describes the 
evolution of Q in the "background" of the radial component, which, in turn, is given by 
the first equation in each case. The evolution of r is easy to understand just looking at 
its equation. First, it falls down as some power of t, which depends on n, and then at 
H Ho starts to oscillate. Numerical comparison of the full system of equations without 
all these simplifications which we made above shows that aX,t<to this approximation is 
reliable. 



6. Conclusions 

Clearly thermal effects are important in the evolution of AD condensate and modify 
the original mechanism of AD baryogenesis. We find that in the n = 1 case it is difficult 
to reproduce the known value of the baryon/entropy ratio even if taking into account 
some uncertainty due to the choice of the value of the reheating temperature, so that 
this case should probably be considered borderline at best. However, the n — 2,3 cases 
give plausible results for a wide range of the parameters of the potential. We have shown 
that the naive estimate of ns/n^ fails by several orders of magnitude to reproduce the 
observed numerical value. This effect is most clearly seen in the n = 2, 3 cases, when 
the radial component of the AD field evolves in the logarithmic potential. This happens 
due to the behavior of the charge before it freezes, which might seem a bit surprising. 
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Instead of decreasing with Ho it actually slightly grows. We have also shown numerically 
that dependence on a few other parameters of the potential such as the phase difference 
between the A terms and the coefficient of the A term b is of the expected form. 
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7. Figures and Captions 



Fig.l: This graph illustrates the behavior of the asymptotic value of 
HBt^ as a function of ln(y) for n = 1. The choice of parameters here is 
a = b = 1; 6b = tt/S. 

Fig. 2: This graph illustrates the behavior of the asymptotic value of 
nst'^ as a function of ln(a) for n = 2. The choice of parameters here is 
a = b = 1; 6b = ir/S. 

Fig. 3: This graph illustrates the behavior of the asymptotic value of 
n^t^ as a function of \n{a) for n = 3. The choice of parameters here is 
a = b = 1; 6b = ir/S. 

Fig. 4 This graph illustrates generic behavior of n^t^ as a function of 
time in the n = 1 case; cases a, b, c, d correspond to ya = 0.01, yb = 0.005, 
yc = 0.0025, yd = 0.001. The values of other parameters are taken to be 
6b = 7r/3, a = 6 = 1 = A. 

Fig. 5: This graph illustrates generic behavior of UBt^ as a function of 
time in the n = 2 case; cases a, b, c correspond to aa = 0.1, ab = 0.07, 
ac = 0.03. The values of other parameters are taken to be 6b = 7r/S,a = 
6 = 1 = A. 

Fig. 6: This graph illustrates the behavior of n^t^ for n = 3, case a 
corresponds to a = 0.1, case b to a = 0.7, and case c to a = 0.5. The 
values of other parameters are taken to be (5?, = 7r/3, a = 6 = 1 = A. 

Fig. 7: This graph illustrates the behavior of the asymptotic value of 
nst^ as a function of sin[6b) for n = 2. The choice of parameters here is 
a = 6 = 1; q;=0.05. 

Fig. 8: This graph illustrates the behavior of charge as a function of b 
for n = 2. The choice of parameters here is a = 1, = 7r/3, q;=0.03 
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Figure 1: 
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Figure 2: 
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Figure 3: 



16 




Figure 4: 
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Figure 5: 
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Figure 6: 
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Figure 7: 



20 




Figure 8: 
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